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Abstract. GX 3+1 is a low- mass X-ray binary that is persistently bright since its discovery in 1964. It was found 
to be an X-ray burster twenty years ago proving that the compact object in this system is a neutron star. The 
burst rate is so low that only 18 bursts were reported prior to 1996. The Wide Field Cameras on BeppoSAX have, 
through a dedicated monitoring program on the Galactic center region, increased the number of X-ray bursts 
from GX 3-1-1 by 61. Since GX 3-1-1 exhibits a slow (order of years) modulation in the persistent flux of about 
50%, these observations opens up the unique possibility to study burst properties as a function of mass accretion 
rate for very low burst rates. This is the first time that bursts are detected from GX 3-1-1 in the high state. From 
the analysis we learn that all bursts are short with e-folding decay times smaller than 10 s. Therefore, all bursts 
are due to unstable helium burning. Furthermore, the burst rate drops sixfold in a fairly narrow range of 2-20 
keV flux; we discuss possible origins for this. 
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' 1. Introduction 

GX 3+1 is one of the first discovered cosmic X-ray sources. 
[ It was det ected during a n ^ero6ee-rocket flight on June 
' 16, 1964 l)Bowver et all fl965b Ever since, GX 3+1 is 
' one of the brightest persistent sources with an average 
2 — 10 keV flux half that of the Crab. Its X-ray intensity 
_ shows a mild variability on an hourly time scale and pos- 
■ sibly a sinusoidal-like v ariation on a time scale of years 
; (|Makishima et al.lll983l) . 

Althou gh an optical c ounterpart has not yet been iden- 
tified (e.g- JNavlor et alJ [l991). GX 3+1 is presumably a 
low mass X-ray binary (LMXB) where a neutron star (NS) 
is accompanied by a low-mass star of spectral type A or 
later. In such a system (for a review, see IWhite et al.ll995l) 
the companion overflows its Roche lobe and dumps matter 
onto the NS through an accretion disc. The accumulation 
of matter on the NS surface leads to pressures and tem- 
peratures that are high enough to initiate nuclear fusion 
processes of the accreted hydrogen a nd helium to highcr- 
Z nuclei ijHansen fc van Hornlll97fj) . The circumstances 
of this process are such that the fusion may be unsta- 
ble; only at accretion rates in excess of the Eddington 
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limit is the burning thought to be completely stable 
(for solar abundances of the accreted matter; Fujimoto 
et al. 1981). Unstable fusion occurs through thermonu- 
clear flashe s which may be observable in the form of X- 
ray bursts ijWooslev fc TaanJll976t iMaraschi fc Cavalierd 
^77) . Such so-called type-I x-ray bursts are characterized 
by a fast rise (one to ten seconds), exponential- like decay 
(e-folding decay times between a few seconds and tens of 
minutes), a black body spectrum of kT ~l -3 keV, and 
coolin g during the decay (for a review, see iLewin et alJ 
For the fla sh to be ob servable, helium is essential as 
a fuel ingredient l)Josglll97'^ . The amount of hydrogen in- 
volved in the b urning layer determines t he longevity of the 
explosion (e.g.- lTaam fc Pickluinlll979|) . Curre ntly, about 
half o f all LMXBs have been seen to burst (e.g.. lin 't Zandl 
l200l|) : GX 3+1 is one of them. 

GX 3+1 was first seen to burst in 1983 (Makishima et 
al. 1983), abo ut eight years after the d iscovery of the first 
X-ray burster ijCrindlav fc Heisell975(l . testifying to an in- 
frequent bursting behavior. Up to now 20 bursts have been 
reported in the literatur e from GX 3+1: 15 with Hakuch o 
llMakishima et alJll983l), 2 with Gin qa l)Asai et alJll993|) . 
1 with Granat (^Pavhns kv et alJl994l) and 2 with the Rossi 
X-ray Timing Explor er (RXTEl ijKuulkers fc VcLn der KlifJ 
l200nl:lKuulkail2002^ . Nineteen of these bursts are short 
with decay times less than 10 seconds. The remaining 
burst, seen with the All Sky Monitor (ASM) on RXTE, is 
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a sup erburst with a decay time of ~1.6 hours ijKuulkerj 
|2002|) . One of the short bursts was seen to exhibit a quick 
(i.e., less than 2 seconds) radius expansion (factor of two) 
phase, indicating that the burst luminosity was near or 
at the Eddington luminosity, causing the NS atmosphere 
to expand due to radiation pressure. This implies a dis- 
tance to the source of 4.5 kpc, for a hydrogen-rich NS 
atmosphere l|Kuulkers fc van der Khsl Eooni. 

In a simple picture a larger accretion rate would imply 
that ignition conditions are reached in quicker succession 
and burst rates should correspondingly be higher, under 
the assumption that the amount of fuel burned per burst 
is similar. Remarkably, the opposite is usually observed. In 
the majority of the roughly ten cases where this issue could 
be investigated (including GX 3-1-1) the burst frequency 
was found to be smaller at higher accretion rates, assum- 
ing that higher fluxes reflect higher accretion rates. Only 
in one case is there clear ev idence for a positive correlation: 
GS 1826-24 l|Cocchi et al.lii20nm . The anti-correlation be- 
tween flux and burst frequency prompted the fuel 'leakage' 
hypothesis: fuel is burnt in a non-obs ervable way with a 
rate t hat increas es with accretion rate l|van Paradiis et aD 
Il988|) . Recentlv. iBildstenI l|200(]() proposed as an alterna- 
tive that the explosion may be occurring on only part 
of the NS surface, and that the expected proportionality 
for the burst frequency is with the accretion rat e per unit 
area instead of the global accretion rate (see also lMarshaJ 
Il982l) . 

We here discuss observations of GX 3+1 that were 
obtained with the Wide Field Cameras on BeppoSAX be- 
tween 1996 and 2002. This data set is unique because of 
its sheer volume, giving rise to the detection of 61 bursts, 
which allows for a more detailed study of the infrequent 
bursting behavior of GX 3-1-1 than was hitherto possible. 
The merits of such a study are increased by the slow but 
large change in brightness over the course of these obser- 
vations, enabling an analysis of the dependency of burst 
properties, such as recurrence times, on flux and possibly 
accretion rate. 



2. Observations 

2.1. BeppoSAX observations 

The Italian-Dutch X-ray satellite BeppoSAX l|Boella et al.l 
Il997,) . launched on April 30, 1996, and switched off exactly 
six years later, c arried two Wide Field Cameras (WFCs; 
l.lager et al.ll997l) . They were active between 2 and 28 keV 
and had field of views (FOVs) of 40° x 40° with an angu- 
lar resolution of 5', and a source location accuracy of 0.7' 
(99% confidence). One of the core programs of BeppoSAX 
was a Galactic center monitoring program with the WFCs. 
The observations were made during two visibility windows 
per year, from mid February to mid April and from mid 
August to mid October (see Table ^ and Fig. . Since 
GX 3+1 is located close to the Galactic center, it was usu- 
ally near the center of the FOV, so that the sensitivity was 
optimum. The complete dataset consists of 12 Galactic 
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Fig. 1. ASM 1.5-12 keV light curve (one-day averages) 
from August 1996 to April 2002. For clarity, data points 
with error bars larger than 1.25 ctss~^ were not plotted 
(i.e., 16% of all points). Below the light curve the time 
spans of the twelve BeppoSAX Galactic center campaigns 
are shown with horizontal bars. Below the time spans the 
burst times arc indicated with vertical lines, and we indi- 
cate the number of bursts per campaign. 



center campaigns, including serendipitous secondary ob- 
servations near the Galactic center, with a total net expo- 
sure time of 6.2 Ms. In this dataset 61 bursts were identi- 
fied with GX 3+1. The bursts were found through analyz- 
ing the full bandpass (2-28 keV) light curves of GX 3+1 
at a time resolution of 5 seconds. The peak fiux detection 
threshold for bursts that have an e-folding decay time of 
at least 3 seconds is then expected to be Ictscm^^s^^, 
or roughly half the fiux of the Crab. 

2.2. ASM observations 

The All Sky Monitor on board RXTE l|Levine et al.lll996() 

monitored GX 3+1 over the six years covering the WFC 
era (and beyond). Fig. ^ shows the 1.5 to 12 keV light 
curve. Plotted are one-day averages of the 90 second dwells 
taken by the ASM. It shows that the intensity of GX 3+1 
slowly changed by a factor of two. At first the source was 
in a 'high state' with a flux of about 0.4 times that of 
the Crab (75 cts s~^). Then the intensity dropped over the 
course of about 20 months until it reached a 'low state' 
with a fiux that is 50% lower. GX 3+1 stayed in this low 
state for about 19 months after which the intensity slowly 
increased until it reached the high state after 10 months 
again. The reason for using the ASM data set is that it 
provides a more uniform coverage than the WFC data set. 
Thus, one can get a better idea how the source behaved 
between bursts when the recurrence time is longer than 
the typical WFC observation. 
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Table 1. Summary of the twelve Galactic center observation campaigns of the WFCs. For every campaign it shows 
the campaign number, year, begin and end time (ibogin and tend), the number of observations, the net exposure time 
(tcxp), the number of bursts seen from GX 3+1 during that campaign, the derived burst rate and the average persistent 
ASM intensity (F^). 
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3. Burst analysis 

A time-resolved spectral analysis was performed of all 61 
bursts. The time resolution is limited by the statistical 
quality of the data, which is determined by the burst flux, 
duration and the off-axis angle. We required a resolution 
so that the 1 sigma errors in the spectral parameters re- 
mained within 20%. 

The 2-28 keV light curves, at a resolution of one sec- 
ond, served as a guide for setting up the time intervals 
for which burst spectra were extracted. The first interval 
was set to start when the photon flux showed a signifi- 
cant trend towards a clear peak. In this regard it is useful 
to note that all 61 bursts showed onsets that lasted less 
than three seconds. The end of the burst was defined as 
the instance when an exponential fit function to the decay 
dropped to the average persistent level within one sigma. 
We also included two time intervals covering the persistent 
emission, namely a 10-15 seconds interval just before the 
burst onset and a 80-90 seconds interval from the end of 
the burst. In total three to five time bins were defined for 
each burst. For each time interval, dead-time corrected 
background-free spectra were extracted in the 2-28 keV 
band. During the analysis we ignored channels 1, 2 and 
31, and included a systematic error of 1.5% per channel, 
to accommodate calibration uncertainties. 

For each burst, we assumed that the spectrum of the 
persistent emission was constant and modeled it with 
thermal bremsstrahlung. This is a simplification with re- 
spect to the more sophisticated models discussed in sec- 
tion 13 but is allowed because 1) the statistical qual- 
ity of the short (115 sec) segments is insufficient to 
discriminate between models, and 2) we are not inter- 
ested in the details of this spectrum at this point but 
merely in (implicitly) subtracting the persistent emission. 
The spectrum of the burst emission was modeled by a 
single-temperature black body, whereby the temperature 
and normalization were left free over all burst time in- 
tervals. On all spectral models, a fixed absorption was 



applied, following iMorrisq n ^ Mc^C^^muu^il 1 198,1 ). with 
TVh = 1.59 X 10^^ cm-^ l lOosterbroek et al.ll2001f) . This 



model fits the data satisfactorily (Xr ranges up to 1.5 with 
at least 80 dof for at least 3 spectra). For each burst in- 
terval the fit provides an estimate for the radius (R) of 
the emission area (i.e., of a sphere and assuming isotrop- 
ical emission), the black body temperature (fcT) and the 
predicted unabsorbed bolometric flux (F) . An example of 
the results of the spectral analysis is shown in Fig.|2| Since 
the time bins are broad, the parameters are not very time 
sensitive. For instance, there is not a high enough tempo- 
ral resohatim_to_detect_a^ such as seen 
bv lKuulkers fc van der KlisI ton(t) . which lasted less than 
two seconds. Also the peak bolometric flux is lower than 
the true peak flux for the same reason. Therefore, in the 
remainder we will consider only the observed peak photon 
flux (-Fpcak), since in this paper we only compare bursts 
seen by the WFC. The total fluence (£'b) is determined by 
integrating the mean observed fluxes over the time bins. 

Fig.Oshows the burst parameters for all 61 bursts plot- 
ted versus time. To search for variability in each parame- 
ter, Xi was determined with respect to the weighted mean. 
This shows that the emission area is consistent with being 
constant (x^ — 0.87). The weighted mean of its spherical 
radius is 4.46 ± 0.08 km for a distance of 4.5 kpc, with 
a standard deviation of 0.44 km. The non- variability sug- 
gests that the emitting radius area is the complete surface 
of the NS and not a variable fraction. However, t he radius 
is sm all with respect to the expected value (e.g.. iHaenseJ 
l20nih . Some of this may be explained by inverse Compto n 
scattering in a hot electron cloud (see lEbisuzakilir987f) . 
Also, we suspect that the distance estimate of 4.5 kpc may 
be too low (see discussion). The other parameters do show 
variability, but only the variation in peak flux is (anti) 
correlated with the state of the persistent flux: when the 
source is in the low state it shows a larger dynamical range. 
The distribution of the peak fluxes in the low state is more 
or less Gaussian with a mean of 2.85 cts cm~^ s~^ and a 
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Fig. 2. Example of the results of the spectral analysis of 
a burst. In the top panel the 2-28 keV light curve of the 
burst is shown. The horizontal dotted line is the deter- 
mined persistent flux. The vertical dotted lines show the 
bin boundaries. The second panel shows the radii of the 
emitting region on the NS (assuming isotropic emission 
and a distance of 4.5 kpc). The radii of the three bins are 
constant over time. The third panel shows the black body 
temperature of the NS. In the bottom panel the average 
bolometric black body flux is shown. 

variance of 0.35ctscm^^ s^^. All the bursts in the high 
state have peak fluxes less than 3ctscm^^s~^. The low 
state distribution predicts a chance probability of 97.5% 
for seeing at least one burst with a peak flux higher than 
3ctscm~^s~^ in the high state. Thus, the chance proba- 
bility of not seeing a burst with a peak flux higher than 
Sctscm^^s^^ in the high state is only 2.5%. Therefore, 
bursts during the low state convincingly reach higher (i.e., 
up to two times) peak fluxes than during the high state. 
We can not say anything conclusive about the decay times 
(r) other than that for all bursts the decay times are 
short (i.e., less than 8 seconds) with a weighted mean of 
3.63 ± 0.10 s. This is similar to the bursts from GX 3+1 
seen with other instruments except for the superburst. As 
mentioned in the introduction, the longevity of a burst 
is indicative of the amount of hydrogen mixture in the 
helium fuel for the flash. The above decay times suggest 
that we are dealing with pure helium burning. In order 
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Fig. 3. In the top panel the ASM 1-12 keV hght curve 
is shown together with the WFC coverage as horizontal 
bars. In the lower panels the burst properties versus time. 
From 2nd to lowest panel: radius of the emission area (as- 
suming 4.5 kpc distance), black body temperature, burst 
fluence, peak flux and decay time. The dotted lines are the 
weighted means of the properties. The radii per burst are 
the weighted means of all time bins, wile the temperature 
is derived from the first time bin since the burst onset. 



to double check this, we averaged all bursts separately for 
the low and the high state, see Fig. 0] Thus, the dynamic 
range of the flux measurements is increased by at least a 
factor of 4. This figure shows that 1) there is no trace of 
prolonged burst emission and 2) the profiles are identical 
within the statistics provided by the data. For compari- 
son, we added in Fig.0]the average WFC-measured profile 
of all long bursts from GS 1826-24. This source is a well- 
docu mented exhibito r of mixed hydrogen/helium flashes 
(e.g., lBildsten ''2000V These bursts are clearly a factor of 
roughly six longer than those of GX 3+1. We conclude 
that all bursts from GX 3+1 are short bursts which are 
thought to be the product of a thermonuclear runaway 
process that is fueled by pure helium. 
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Fig. 4. Comparison of average 2-28 keV burst profiles. 
From top to bottom, the average profiles of GX 3+1 in the 
low state, GX 3+1 in the high state and GS1826-24 are 
given. The average decay times are 4.23±0.13 (x? = 1-15), 
4.97 ± 0.28 (x? = 0.61) and 29.21 ± 0.97 (x? = 1.13) 
seconds, respectively. 



4. Burst rate 

We have investigated the burst rate as a function of the 
average persistent intensity (-Fpers)- To assess the expo- 
sure times during which the sensitivity was sufficient to 
detect bursts with peak fiuxes in excess of 1 ctss~^ cm~^ 
(or 0.5 Crab) we excluded all the observations whenever 
the source was located within 30 pixels from the edge 
of the FOV. This amounts up to 6% of the total expo- 
sure. For the remaining data we derived the net exposure 
times, taking into account Earth occultations and pas- 
sages over the South Atlantic Geomagnetic Anomaly (see 
texp in Table^l. The burst rate was calculated by dividing 
the number of bursts by the net exposure time for each 
campaign. The 68% confidence margins were determined 
by assuming Poisson statistics. This results in asymmet- 
ric errors. As no bursts were seen during the third and 
the twelfth campaign only upper limits to the burst rates 
could be derived. They are 3.8 and 12.5 Ms'^ (68% con- 
fidence), respectively. For the complete list we refer to 
Tabled Since the persistent emission and the burst rate 
of the different campaigns are sometimes similar, we av- 
eraged them over several campaigns. Campaigns 5, 6 and 
7 were averaged for the low state and campaigns 1, 2, 
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Fig. 5. Burst rate plotted versus the 1.5-12 keV ASM 
intensity. The change in burst rate is a factor of six. The 
scale on the right hand side is the equivalent recurrence 
time for the specific burst rate (d for days and h for hours). 
The dashed lines are the best fit of the high- and low 
burst rate level. The dotted line is the best straight line 
fit. In the upper right corner relevant information is given, 
namely the number of bursts, the net observation time 
and the number of campaigns for each data point from 
left to right. 

9, 10, 11 and 12 (a null observation) for the high state. 
Campaigns 3, 4 and 8 are the campaigns in which the over- 
all intensity decreased or increased, of which campaign 3 is 
a null observation. For every campaign we calculated the 
average persistent intensity by averaging the ASM mea- 
surements over the complete time span of the campaign. 
Thus, biases due to the variability on time scales shorter 
than the burst recurrence time were suppressed. In Fig. [3 
the burst rate is plotted versus the persistent emission. 
At first sight, one can discriminate two levels of burst 
rate, with a quick transition from a high burst rate to 
a low burst rate at an intensity between 18.3 cts s~^ and 
20.5 cts s~^. To test models for the behavior of the burst 
rate as a function of intensity one may not use the stan- 
dard x^-minimization, because the statistical fluctuations 
are not expected to be distributed as a Gaussian. Instead 
we chose to use the statistic 



N 



{di -rriiY 



where di is the i-th measurement, rrii the pre- 
dicted measurement and N the number of measurements 
{N = b). S was minimized through a grid search. Chance 
probabilities were evaluated through Monte Carlo simu- 
lations, by counting the number of times that simulated 
S'-values were larger than measured S'-values. 

To find the most probable behavior three models were 
investigated: a constant, a linear and a step function. For 
the constant function no satisfactorily fit could be found, 
with a chance probability of 9 x 10^^. For the linear model 
a better fit was found. For a function of the form a + bx 
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Fig. 6. Burst rate plotted versus the 2-28 keV WFC flux. 
In essence is this the same plot as FigEl but now the per- 
sistent flux is averaged over smaller time intervals. Every 
campaign has been plotted, because the spread in persis- 
tent intensity is larger. One sees that the transition be- 
tween the high and low burst rate may be broader than 
in Fig|Sl(i.e., 20% in stead of 10%). 



the values a = 46.46 Ms~^ and b = -1.69 Ms~ Vets s"^ 
were found. The chance probability is 21.6%. The best fit 
is for the step function. The two levels that were found are 
21.1±3.0Ms"i and 3.67±1.0Ms-^ The chance probabil- 
ity is 95.3%. Fig.|Slshows the two fits. If the true relation- 
ship is a linear function, the chance that in a measurement 
like ours the minimum S is smaller for a fit with a step 
function than for one with a linear function is 7.2 ± 0.2%, 
as shown by Monte Carlo simulations. This is the prob- 
ability that we are dealing with a linear function rather 
than a step function. 

In order to obtain an idea how the transition depends 
on the definition of the averaging of the persistent flux, 
we plot in Fig. the burst rate as a function of the WFC 
average persistent flux. It shows that the transition be- 
tween the high and low burst rate is less narrow. Also the 
distribution of the fluxes during the high state is broader. 
This is the reason why we have not plotted averages over 
campaigns like in Fig. |S1 The broader transition is also 
apparent if we plot versus ASM flux if averaged over the 
times when the WFCs were on target rather than the com- 
plete time span of each campaign. 

5. Persistent emission and the mass accretion rate 

At time scales larger th an one minute, GX 3-1-1 exhibits 
two kinds of variability IjAsai et al One has a time 

scale of order minutes to hours and the modulation depth 
in the 1-20 keV intensity is typically 20%. The inten- 
sity either correlates with the spectral hardness (when the 
source is said to be in the 'banana' state) or not ('island' 
state). The 7-18/1-7 keV hardness ratio is observed to 
change by about 50%. The other kind has a time scale of 
order years and is what causes the source to exhibit the 



so-called 'low' and 'high' states (Fig. P). In the following 
we will concentrate on the latter (slow) variability because 
the minimum wait time between bursts is longer than the 
time scale of the faster variability; also the faster variabil- 
ity is presumed to be unrelated to changes in the mass 
accretion rate (Asai et al. 1993). The slow variability has 
been shown nicely by 19 87-1991 measure ments with the 
all sky monitor on Ginga ijAsai et al.lll993l) . The trends in 
those observations are quite similar to the ASM measure- 
ments. The same change in intensity is observed (a factor 
of 2); such a variation has bee n observed in earlier times as 
well (|Makishima et alj[l98^ . The Ginga, ASM and ear- 
lier data sets even suggest that the flux oscillates semi- 
sinusoidally with a period of 6 to 7 years. Interestingly, the 
Ginga measurements with its Large Area Detector show 
a const ant 1-20 keV spe ctral shape over the low and high 
states l|Asai et al.lll993|) . 

The stability of the X-ray spectrum over the low and 
high states is more or less confirmed by the WFC mea- 
surements in the 2-28 keV band in the sense that spec- 
tral changes are minimal, see Table |21 and Fig. [3 We 
made a selection of one long observation (i.e., longer than 
one day, except for the last campaign which did not in- 
clude such a long observation) per campaign and inves- 
tigated the full-resolution spectra. A general fit with an 
absorbed power law plus disk black body, fix ing Nn again 
at 1.59x10 cm~ (lOosterbroek et al .l200lh and coupling 
the disk black-body temperature and power law index over 
all 12 spectra, yields x^/'^ — 1-31 {i^ ~ 310; assuming a 
systematic uncertainty of 1.5% and ignoring the insuffi- 
ciently calibrated spectral channels 1, 2 and 31 of each 
spectrum). If the temperature of the disk black body is 
allowed to vary over the 12 spectra, the fit improves to 
X^/i^ — 1-027 {v — 299). Uncoupling the power law in- 
dex instead of the temperature yields a slightly worse fit 
with Y^/i/ = 1.06 (f = 2 99). Applying the same model as 
lOosterbroek et alJ ^Mh) on the 0.5-200 keV NFI spec- 
trum, leaving free only the fluxes of the 3 components, the 
optical depth and the black body temperature, but keep- 
ing them coupled over the 12 spectra except for the disk 
black body normalization, results in a reasonable though 
formally not acceptable fit: x^/'^ = 1.354 {v = 320). When 
the disk black body temperature is allowed to vary, the 
fit improves to x^/^ = 0.945 {i^ = 309). Irrespective of 
the two models, the black body temperature changes only 
modestly by about 15% over the 12 spectra. 

We conclude that the one-day averaged spectral shape 
shows only a very modest change over a time scale of years. 
This change is completely attributable to the disk black 
body component. This black body component is likely due 
to emission from the inner edge of the accretion disk and 
is pr obably driven by changes in the mass accretion rate 
re.g..lRa,rretjE00T]). 

A parameter that is more meaningful than the 2-12 
keV photon count rate as measured with the ASM is 
the unabsorbed bolometric flux fboi because that should 
approximate the mass accretion rate more accurately. If 
the spectral model is extrapolated to beyond the WFC 
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Table 2. Results of best-fit spectral models for selection of parameters. The fits were applied to twelve WFC spectra 
of GX 3-1-1, one per campaign. 



W (2) (3) (4) (5) (6) (7) (8) 



■iDSI Q 77- 


-50'H90 86 


10.7 


9.14 


2 


291+0 


013 


91 0+0 6 


9 1 74+0 093 


21.0+1.0 


13.7 




-505^19 OS 


10.6 


9.07 


2 


255+0 


010 


93 1 +0 6 


9 1 40+0 093 

jLi . 1 '-+1 } 1 \ 1. KJ^rJ 


22.1+1.1 


13.6 


^n79'i 4-4-- 


-50796 54 


8.95 


7.74 


2 


234+0 


012 


90 0+0 5 


9 081 +0 097 


19.9+1.2 


11.6 


50910.10- 


-50911.25 


6.31 


5.46 


2 


155+0 


019 


15.2+0.6 


1.861+0.041 


18.2+1.8 


8.3 


51057.24- 


-51058.38 


6.03 


5.15 


2 


127+0 


026 


14.2+0.8 


1.768+0.044 


20.3+2.2 


7.8 


51270.08- 


-51271.21 


6.60 


5.62 


2 


211+0 


020 


13.9+0.6 


1.938+0.041 


16.2+1.5 


8.5 


51422.68- 


-51423.81 


5.95 


5.19 


2 


109+0 


020 


15.9+0.8 


1.779+0.045 


20.1+2.2 


7.9 


51618.05- 


-51619.21 


8.33 


7.22 


2 


158+0 


016 


20.4+0.7 


1.966+0.029 


22.6+1.5 


10.9 


51802.28- 


-51803.42 


7.98 


7.87 


2 


256+0 


Oil 


17.2+0.4 


2.089+0.032 


16.4+1.1 


10.3 


51989.20- 


-51989.88 


11.5 


9.76 


2 


336+0 


022 


21.2+1.1 


2.240+0.030 


20.3+1.2 


14.6 


52181.79- 


-52182.83 


11.1 


9.54 


2 


299+0 


010 


22.8+0.5 


2.203+0.023 


21.0+1.0 


14.3 


52353.24- 


-52353.71 


9.72 


8.34 


2 


293+0 


016 


20.1+0.8 


2.179+0.029 


18.3+1.1 


12.5 



(1) Observation time span (MJD); (2) unabsorbed 2-28 lieV flux (in units of 10 " erg cm ^s (3) unabsorbed 2-10 IceV 
flux (in units of 10"^ erg cm"^s"^); (4) diskbb ItTin (keV); (5) diskbb Rf^cos(e) (km^); (6) and (7) as (4) and (5) except in 
combination with a Comptonized spectrum instead of a power law; (8) unabsorbed 0.01-100 keV flux predicted from the latter 
model (in units of 10~^ erg cm~^s~^); 



energy is released in another form than radiation. This 
emphasizes that our lower limit for the accretion rate is a 
firm one. We note that the mass accretion per unit area 
m can be estimated ind ependently from the distance by 
{FboiR/GM){d/ R^ppf) l)Bildstenll2OO0l) . For the neutron 



star canonical values 



is above 2 x 10^ g cm ^s ^. 



Channel energy ^keV) 



Fig. 7. Examples of incident photon spectra during low 
state and high state, as measured with the WFC on MJD 
51057.24-51058.38 and MJD 52181.79-52182.83 respec- 
tively. 



bandpass, Fboi is estimated to vary between 8x10^^ and 
14x10^^ erg cm~^s~^ (see Table 121. Such an extrapola- 
tion is j ustified because a s t udy o f the 0.5-200 kcV spec- 
trum bv lOosterbroek et alJ l|200lh . taken on MJD 51420- 
21, showed that 99% of the 1-200 keV flux is contained 
within 1-20 keV. For a distance of 4.5 kpc the flux trans- 



lates to a bolometric luminosity of 2-3x10'^^ erg s^^ (7- 
12 X lO'^^ erg s~^ if at a canonical distance of 8.5 kpc, equal 
to that of the galactic center). The smallest observed flux 
(8x10^^ erg cm~^s~^) in combination with the smallest 
possible distance (4.5 kpc) and the largest possible NS 
radius (10 km), while assuming a NS mass of 1.4 M0 
and a luminosity given by G M M/R, yields the minimum 
mass accretion rate in GX 3+1: 1.7x10^^ Mq yr^^ (or 
10% of Eddington for a hydrogen-dominated NS atmo- 
sphere) . It is likely that part of the liberated gravitational 



6. Discussion 

Our observations of the burst properties as a function of 
persistent flux can be summarized as follows. Within a 
narrow range of persistent flux of about 10% (i.e., 18.3 
to 20.5 ASM ctss^^), the average burst frequency drops 
by a factor of 6 from once per 13.2 hours to once per 
3.2 days. Though the best fit is a step function, we can 
not rule out a gradual transition. This is the first time 
that GX 3+1 is seen to be burst active in the high inten - 
sity state (c.f.. iMakishima et alJll983l: lAsai et al.lll993|) . 
From all other burst properties only the peak flux seems 
to be (anti) correlated with the persistent flux. It occa- 
sionally shows higher values, by a factor of up to 2, dur- 
ing the low state. In terms of the decay time, there is no 
difference between the low and high state (i.e., they are 
short) and all bursts appear to be predominantly fueled 
by helium. As noted in the introduction, an anti corre- 
lation between persistent flux and burst frequency is not 
unexpected. However, the issues that set GX 3+1 apart 
are the quickness of the transition and the fact that on 
both sides of the transition only helium bursts are seen. 
There is perhaps one other source, that we know of, that 
may show a similar effect but the num ber of bursts seen i n 
the high state is only two (4U 1636-53: lLewin et al.lll98'^ . 
The origin of the anti correlation is not well understood. 
The fact that in GX 3+1 the anti correlation takes place 
over a rather small range of persistent flux may perhaps 
be an important clue. 

The discrete-like nature of the transition in burst rate 
is reminiscent of the boundaries predicted by theory be- 
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tween re gimes of unstable n uclear burning on accret- 
ing NSs ijFuiimoto et al.l[l98l|) . For mass accretion rates 
M > 2 X 10""'^° MQyr~^, the accreted hydrogen-rich ma- 
terial is expected to undergo stable hydrogen burning. At 
the bottom, a helium layer builds up in which the pressure 
increases because of the continuing accretion. At a certain 
moment the density becomes high enough to ignite the 
helium in an unstable burning process giving rise to a 
thermonuclear flash. If M < (4 - 11) x lO"^°M0yr"\ 
the pressure is high enough for helium ignition only in 
the helium layer; however for higher M's, the bottom of 
the hydrogen layer may also have sufficient pressure. As 
a result, the flash takes place in a mixture of helium and 
hydrogen and proton captures occur by the nuclear prod- 
ucts of the helium burning which makes the flash extended 
because of the longer ti me scale of the resulting (3 decays 
l)Taam fc Picklumll979(l . Below 2 x 10^^° Mq yr-\ the ac- 
creted hydrogen no longer undergoes stable burning, but 
unstable burning may occur at lower temperatures be- 
cause of pressure build-up by the accumulating matter. 
The heating produced by this runaway process may ig- 
nite the helium accreted from the companion, provided the 
mass of the accumulated layer is large enough. In general, 
the amount of hydrogen participating in the flash is larger 
in this regime than in the regime Af > 2 x 10~^° M© yr~^. 
Finally, there is a regime above roughly lO^^Moyr"^ 
where even the helium is burned in a stable manner. This 
is near the Eddington limit and may not be observable. 

The above regimes can experimentally be confirmed 
through measurements of burst durations. All observable 
bursts involve helium burning. For pure helium, the burst 
durations are predicted to be of order 10 sec, while the 
rise times are shorter than two seconds. Whenever hy- 
drogen participate s , the burst duration becomes longer. 
Ivan Paradiis et all ^|l988^ confirmed the correspondence 
between theory and measurement, although for a limited 
data set. A particularly nice i llustration of the regira es in 
a single source is 4U 1705-44 l|Langmeier et al.lll987|) . 

In GX 3-1-1, the minimum M that is consistent with 
the persistent flux is 1.7 x 10~^ M0yr~^. This is near 
or slightly above the upper boundary of the range where 
pure helium bursts are expected to occur. Therefore, the 
data appear discrepant with theory, even if abundances 
are sub-solar. Still, the factor-of-two dynamic range of 
the persistent flux in GX 3-1-1 fits comfortably in the 
factor-of-five dynamic range predicted by theory for he- 
lium bursts. How do other bursters behave that have 
similar or higher persistent luminosities? Only a minor- 
ity of bright LMXBs burst; we can think of the fol- 
lowing four cases. Cyg X-2 is a Z-source with a lu- 
minosity close t o Eddin gton. Thus far, 18 bursts were 
detected llKahn fc Grindla^ Il984l: iKuulkers et all Il995t 
ISmalfl Il998l iTitarchuk fc Sha,noshnikovl l2002l) . All are 
short, like in GX 3-1-1. The onl y other Z-sou rce that 
exhibits bursts is GX 17+2 Ce.g.. iKuulkers et al. ,2002) . 
In that case seven out of ten bursts are clear examples 
involving hy drogen burning. Other sourc es are Cir X-1 
i Tennant et al. .1 98fia.M and GX 13-^1 l|Matsuba et a,lJ 



While Cir X-1 shows somewhat longer bursts, the 
two bursts detected from GX 13-1-1 are again short. We 
propose that the regime of mixed hydrogen/helium bursts 
at high accretion rates is seldom observed, that the regime 
with pure helium bursts may extend up to higher mass 
accretion rates, and that the transition of burst rate in 
GX 3-1-1 is not due to change of regime. 

The usual explanation for the anti-correlation is that 
a larger part of the helium is burned in a stable manner in 
the h igh state i Makishima et alJll983Hvan Paradiis et alJ 
Il988|) . Thus, less fuel is burnt in thermonuclear flashes. 
This is supporte d by the measure ment of the burst pa- 
rameter a (e.g., iLewin et al.1ll993|) which is the ratio of 
the fluence of the persistent emission from burst to burst 
and that of the burst. For the low state the average value 
is 1.7 X 10^ and for the high state 2.1 x IC*. Both these 
values are much higher than what is expected and is nor- 
mal for most X-ray bursters (between 25 and 250). It is 
a strong indication that stable helium and/or hydrogen 
fusion is happening in b oth the low and the high state 
l)van Paradiis et alTll988l) . From the ratio in a the sug- 
gestion is that the energy production of the stable fusion 
which goes into radiation increases 12-fold from the low 
to the high state. Stable helium burning is predicted to be 
possible only at near-Eddington accretion rates or above, 
while GX 34-1 accretes at most at 30% of Eddington (see 
below). However, this difficulty could be circumvented by 
assumin g that accretion takes place on a con fined area of 
the NS iFushiki fc Lamblll987tlB'ildstenll200f]l) . Anyway, if 
stable helium fusion is the process responsible for the anti- 
correlation, our detection of the rather discrete nature of 
this anti-correlation should provide an important test to 
theory. 

There are the open issues of the higher peak fluxes 
during the low state, the comparatively small distance for 
a source that has such a small angular separation from 
the galactic center, and the doubtful small radii for the 
burst emission areas. We prop ose that they may all b e re- 
solved in the following manner. ' Ebisuzaki et all l|l983|) de- 
termined that radiation from helium bursts may be intense 
enough to expel the hydrogen atmosphere, thereby leav- 
ing a hydrogen-poor photosphere. Such a photosphere has 
a higher Eddington limit so that luminosities may reach 
higher values by roughly a factor of two. This is more likely 
to happen if the hydrogen-dominated atmosphere is thin- 
ner. This is precisely what applies during the low state: 
the accretion rate is twice as small and the bursts recur 
six times quicker so that the accumulated mass from the 
companion between bursts is roughly twelve times smaller. 
Thus, this provides an explanation why peak fluxes are 
higher during the low state. The same explanation has 
recently been applied to the bim odal peak flux d i stribu - 
tion of bursts from GX 354-0 bv iGallowav et alJ \2m± . 
Unfortunately, in our case we are not able to confirm the 
radius-expansion phase in the brighter bursts due to in- 
sufficient statistics. 

However, the only radius-expansi on burst from 
GX 3-hl detected bv lKuulkers fc van der Klis, (200Q1 has a 
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peak flux of 2.3 Crab and happened during the low state. 
The peak flux compares well, within a few tens of per- 
cent, with the maximum peak flux that we find during 
the low state. The above reasoning implies that this burst 
reaches the Eddington limit for a helium-rich photosphere. 
iKuulkers fc van der Klij l|2000tl determined that in that 
case (X=0) the distance would be 6.1 ±0.1 kpc. This, in 
comparison to the 4.5 kpc, is substantially closer toward 
the galactic center (at 8.5 kpc) and seems more plausible. 
Furthermore, the radius of the burst emission scales up 
accordingly to a more plausible value of 6.05 ± 0.11 km 
(weighted mean). The luminosity of the persistent emis- 
sion would increase accordingly. The maximum luminosity 
would be roughly 6 x 10^^ erg s~^, or 30% of Eddington. 
We note t hat 6.05 km is still smaller than the expected 
NS radius llHaense]ll200l . This may be explai ned by in- 
verse C ompton scattering effects as discussed bv lEbisuzakil 
l)l987(l which may introduce a factor of 2 decrease in the 
apparent radius for the black-body temperatures seen in 
GX 3+1. 
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